Plasmons enhance near-field radiative heat transfer for graphene-covered dielectrics 
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It is shown that a graphene layer on top of a dielectric slab can dramatically influence the ability 
of this dielectric for radiative heat exchange. Effect of graphene is related to thermally excited 
plasmons. Frequency of these resonances lies in the terahertz region and can be tuned by varying the 
Fermi level through doping or gating. Heat transfer between two dielectrics covered with graphene 
can be larger than that between best known materials and even much larger at low temperatures. 
Moreover, high heat transfer can be significantly modulated by electrical means that opens up new 
possibilities for very fast manipulations with the heat flux. 

PACS numbers: 44.40.+a, 42.50.Lc, 78.67.-n 
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Radiative heat transfer (RHT) in vacuum at small dis- 
tances between bodies is much increased in the near-field 
regime as compared to that given by the black body law 
[lH3|. It happens due to interaction of evanescent waves 
at distances small in comparison with the thermal wave- 
length Xt — hc/T (here ks = 1)- Particularly strong 
enhancement occurs when bodies can support surface 
modes such as plasmon-polaritons and phonon-polaritons 
[3, [H| • This effect can be used to improve performance of 
near- field photovoltaic devices Q , in nanofabrication 0] , 
and in near- field imaging systems @. 

Graphene attracted recently enormous attention as a 
two dimensional carbon material with unusual electronic 
properties [t| . It is considered as a promising material for 
the development of high-performance electronic devices 
fiol [ill ] . Plasmons in graphene show favorable behavior 
for applications such as large confinement, long propagat- 
ing distances, and high tunability via electrostatic gating 
[12| . In contrast with nobel- metals the plasmon frequen- 
cies lie in the terahertz region that is interesting for ra- 
diative heat transfer, but the topic was not explored yet. 
Heat transfer was considered [13[ only between closely 
spaced graphene and Si02 substrate where plasmons do 
not play significant role. 

Pristine graphene at zero temperature does not sup- 
port plasmon excitations but doped material does [14j . 
On the other hand at finite temperature plasmons exist 
even for undoped material [15| . These thermoplasmons 
were shown to change significantly the thermal Casimir 
force for suspended graphene [l6| and graphene-covered 
materials [17} . In this paper we show that plasmon exci- 
tations in graphene have striking effect on the near-field 
RHT between bodies if at least one of them is covered 
with graphene. 

Usual local materials have fixed frequencies of phonon- 
polariton or plasmon-polariton resonances. Graphene is 
essentially nonlocal material and its plasmon frequency 
changes with the wavenumber. Moreover, it varies sig- 
nificantly with the doping level. These properties make 



plasmons in graphene a convenient tool to control the 
heat transfer between bodies. 

To evaluate the RHT between two bodies 1 and 2 



one has to know the reflection coefficients 



and 



for each body as functions of the frequency uj and the 
wave vector q. These coefficients arc different for each 
polarization fi = s (transverse electric) or /i = p (trans- 
verse magnetic). If the separation d between two paral- 
lel plates is small in comparison with the thermal wave- 
length, d -C Ay, only evanescent waves will contribute to 
the heat transfer. Defining the heat transfer coefficient 
(HTC) as h = $(T,T - AT) /AT for AT -> 0, where 
$(Ti, T 2 ) is the heat flux and AT = T\— T 2 is the tem- 
perature difference between bodies, one finds h = h p + h s , 
where 
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Here the integration variable in the physical terms is x = 
2d v / q 2 ~uj 2 /c 2 . 

Suppose that the body i is a dielectric substrate with 
the dielectric function Si(ui) covered with a layer of 
graphene. Reflection coefficient for this body can be 
found describing graphene as a sheet with a current. 
Then for p-polarization one has [18]: 
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where fco — y/ w 2 /c 2 — q 2 and ki — y SiuS 1 / 'c 2 — q 2 are 
the normal components of the wave vectors in vacuum 
and in the substrate, respectively, and a(q, uj) is the two- 
dimensional (2D) dynamical conductivity of graphene. 
We omitted the polarization index in (J2|) because the 
contribution of graphene can be found using only p- 
polarization. Moreover, similar to the situation with 
the Casimir force [H, [TtJ (see also [l9j]) we can use 
the non-retarded approximation to calculate the effect 



2 



of graphene. Retardation and the contribution of s- 
polarization both are suppressed at least by the fac- 
tor vp/c ~ 1/300, where up is the Fermi velocity in 
graphene. However, it has to be stressed that the non- 
retarded limit can be applied only to the graphene con- 
tribution [17(. For the HTC this contribution is defined 
as 



Ah = h p (r 1 ,r 2 ) - h p (r 10 ,r 20 ), 



(3) 



where h p must be understood as a functional defined by 
Eq. ((T|) and is the reflection coefficient of the body 
i without graphene. To calculate any term h{r\,r 2 ) or 
M r iO: r 2a) separately one has to include both the polar- 
ization and retardation effects. 

Dielectric function of graphene on the interface of 
vacuum and body i can be calculated in the random 
phase approximation Q as e(q,w) = 1 + v c (q)H(q,uj), 
where v c = 2ire 2 /niq is the 2D Coulomb interaction, 
Ki is defined by the environment of the graphene layer, 
2m = £;(0) + 1, and H(q, uj) is the 2D polarizability given 
by the bare bubble diagram. The latter was calculated 
in many papers 0, H3 i here we are using the result [I?} 
that can be applied for both finite temperature and finite 
Fermi level. One can equally express the result via 2D 
susceptibility = (e 2 /q 2 )H(q, w), which, in turn, is 

expressed via 2D conductivity as — iuix(q,w) = 
In this way the reflection coefficient ri in the non-retarded 
limit (c — > oo) can be expressed via the dielectric func- 
tions of the substrate £j(w) and graphene e{q,ui): 



e t - 1 + 2Kj{e - 1) 
Si + 1 + 2 Kl (e - 1) 



(4) 



Random phase approximation describes collisionless 
electron gas, but this approximation is not sufficient 
for the RHT. This is because the graphene contribution 
given by Eqs. (|3]) and ([I]) disappears in the limit Ime — > 
and finite dissipation is of principal importance. At low 
frequencies the main dissipation in _graphene is due to im- 
purities and defects for electrons [12j . It can be included 
via the finite relaxation time r = 7 -1 as proposed by 
Mermin [2l|, where 7 = 10 12 — 10 13 rad/s. In this ap- 
proach the dielectric function s T (q, ui) is expressed via the 
collisionless function e(q,u) taken at complex frequency 
u) — > ui + ij as 



e T (q,ui) = 1 + 



(uj + ij) [e(q, w + ry) — 1] 
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Function e r has to be used in Eq. (f4]) instead of s. 

At finite doping and finite temperature the function 
s(q, uj) can be presented in an analytic form in the limit 
Q = Hv F q/2T < 1 pi}. In this limit the intraband 
transitions dominate in the dielectric function, while the 
interband transitions are suppressed by the factor Q and 
the dielectric function is 



e(q,u) = 1 + 
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where ep = Ep/T is the dimensionless Fermi level and 
a g = e 2 /nhvp is the coupling constant in graphene. 
Function G(6f) is defined as: 



G(e F ) = / dtt 



1 + cosh t cosh e p 
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For small and large values of ep it takes the asymptotic 
values G(ep — > 0) = 2 In 2 and G(ep — > 00) = ep, respec- 
tively. 

Let us assume first that the substrate of the body 1 
has no optical activity in the terahertz range. It means 
that £i(w) ~ £i(0). In this case the reflection coefficient 
has a pole when e(q,uj) = 0(7-) 0). This pole describes 
plasmon in graphene with the dispersion relation 
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When ep — > we reproduce Vafek's thermoplasmon [T5| . 
In the limit of large ep the plasmon frequency does not 
depend anymore on T and we find the plasmon that 
emerging at finite doping 14 1. 

For the best heat transfer the resonances in the op- 
posing bodies have to match each other 4]. When one 
body is covered with graphene one can always find a 
value of the wavenumber q that gives the plasmon reso- 
nance matching the resonance in the opposing body.This 
simple principle gives qualitative explanations for rich 
physics that can be realized between bodies covered with 
graphene. Let us illustrate this statement. 

Suppose that the second body can be described at low 
frequencies by a single Lorentz-Drude oscillator with the 
dielectric function 



E2M = £00 1 + 
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where ui r is the resonance frequency, r <C uj r is the res- 
onance width, A is the amplitude, and £oo is the high- 
frequency dielectric constant. Reflection coefficient r 2 
has the resonance at frequency corresponding to the sur- 
face wave excitation that is determined by the equation 
£2(0;) = — 1. This frequency is 
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Imaginary part of r 2 is shown in Fig. HJa) by the dashed 
line. Plasmon resonances in graphene give peaks in Imri 
that are shown for x = 1 and x = 3. Because x = 2dq 
there is a value of q when frequencies of the resonances 
in both bodies match each other. 

Integration over x in Eq. (TTJ) gives the spectral density 
of the HTC shown in Fig. QJb) for two values of the 
Lorentz resonance. This density is concentrated around 
the resonance in the body 2. Integrating over u> one finds 
the effect of graphene in the HTC. In the rest of the paper 
we present the results for the scaled HTC defined as 



AR = (d/100 rnn) 2 ■ Ah. 



(11) 
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FIG. 1: (Color online) Inactive dielectric covered with 
graphene (body 1) against a dielectric without graphene de- 
scribed by the Lorentz model (body 2). (a) Resonances in 
the reflection coefficients. The dashed line (green) shows the 
fixed resonance in Imr2. Plasmon resonances in the body 1 
are shown for x = 1 and x — 3 together with the enveloping 
line (black), (b) Spectral density for two values of the Lorentz 
resonance, (c) AR as a function of distance. Curves 1, 2, 3 
correspond to ef =0, 5, 10, respectively, (d) AR as a function 
of the Fermi level. 

In Fig. QJc) AR is shown for the Lorentz resonance 
Lj r i = 0.03 eV and three different values of the Fermi 
level. We can compare h — AR ■ (100 nm/d) 2 with the 
black body coefficient w 6 W/m 2 K. One can see that 
graphene provides a significant contribution to the HTC. 
Dependence of AR on the Fermi level is nontrivial and 
strong. It is shown in Fig. (Ud). Value of AR varies 
more than 10 times when the Fermi level changes from 
to 10T. Moreover, this dependence changes from de- 
creasing to increasing when u r changes from 0.03 eV to 
0.1 eV. Indeed, the characteristic plasmon frequency uit* 
is given by Eq. © for q ~ 1/2(1 If ujf > 0J r for all E F 
then AR will decrease when the Fermi level increases. 
In the opposite case AR will increase with Ep while the 
condition oj p < ui r holds true and change to decreasing 
when the condition is broken. 

Effect of graphene is especially strong for identical in- 
active bodies covered with graphene (see Fig. dta)). This 
is because the plasmons frequencies in both bodies coin- 
cide for every q so that the spectral density becomes much 
wider than in the case of one graphene-covered body. 
This density is shown in Fig. [U[b) for different values of 
Ep. The amplitude and width of the curves is essentially 
controlled by the thermal factor in Eq. (fT]). Value of AR 
becomes smaller when the substrate dielectric constant 
£i(0) = £2(0) = So increases. It is also reduced very fast 
when the substrates are not the same. This is because 
the plasmon frequencies do not match any more due to 
different coupling constants a g in the graphene layers of 
different bodies. When one or both of the substrates are 



FIG. 2: (Color online) (a) AT? as a function of ef for inactive 
dielectrics both covered with graphene. Curves 1, 2, 3 are for 
£0 =4, 7, 12. Curve 4 is for ei(0) = 4 and e 2 (0) = 7. (b) 
Spectral density for e = 4 and e F =0(1), 5(2), 10(3). (c) 
Two identical substrates described with the Lorentz model; 
one is covered with graphene; ef =0(1), 5(2), 10(3). (d) Sub- 
strates are the same as in (c) but both covered with graphene. 
Dependence on ef is shown for d = 50(1), 100(2), 150(3) nm. 

metals the effect of graphene disappears, AR — > 0. It 
happens because in the limit £, — > 00 the reflection coef- 
ficients with and without graphene coincide, r$ — > r^o- 

Suppose that body 1 is covered with graphene but both 
substrates are optically active and can be presented by 
the same Lorentz oscillator. In this case the body 1 sup- 
ports two surface waves with the frequencies 
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These frequencies do not coincide with W2s- Only in the 
limit Up — > one has (w is )i — > u)2 S - It means that the 
resonances in different bodies never match each other 
and the term h p (ri,r2) in Eq. ([3]) cannot be large. On 
the other hand, the term /i p (rio, ^20) must be large be- 
cause without graphene on body 1 the resonances coin- 
cide. Therefore, for identical optically active substrates, 
graphene on one of the bodies will reduce the heat trans- 
fer, AR < 0, as demonstrated in Fig. [2jc). Indeed, the 
total HTC given by Eq. (fT]) is always positive. 

When both active substrates are covered with 
graphene then the resonances in different bodies match 
and the spectral density is wide similar to that shown 
in Fig. [2jb) . Difference with the case of inactive sub- 
strates is that the term h p (r 10 , r2o) in Eq. §3§ is nonzero 
and gives important negative contribution. Now AR can 
be negative or positive depending on the parameters as 
shown in Fig. [2jd). 

For applications in nanoelectronics it is interesting to 
have a device that can open and close a heat transfer 
channel at high switching rate [22[ . Bodies covered with 
graphene allow deep modulation of the heat flux at very 
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FIG. 3: (Color online) (a) AT? between two inactive dielectrics 
(eo = 3) covered with graphene as a function of the Fermi level 
difference AEf = Efi — Ef2', 7 is in units 10 12 rad/s. The 
inset shows how the maximal value depends of the distance d. 
(b) Ratio of the HTC for graphene-covered dielectrics (AR) 
and two SiCb plates (T?s;o 2 )-The inset presents dependence 
of the maximal AT? on the relaxation frequency. 

high frequencies. The switching frequency / must be 
smaller than important frequencies in the graphene di- 
electric functions, / -C vf/^ttcI ~ 1 THz but actual re- 
striction comes from graphene electronics that can be 
switched at frequencies up to 100 GHz (23[. The largest 
heat flux one can get for identical substrates inactive in 
THz range that have low dielectric constant in this range. 
Many polymeric materials used in electronics [24[ meet 
these conditions, for example, polyimide. 

When graphene layers on both substrates are equally 
doped, then the HTC is large but does not change much 
with the Fermi level as shown in Fig. [21a). However, 
when the doping in each layer is different the HTC de- 
creases rapidly because of mismatch of the plasmon fre- 
quencies in the graphene layers. Due to this dependence 



one can modulate the heat transfer between bodies by 
changing the difference between the Fermi levels in dif- 
ferent bodies. Situation is illustrated in Fig. 03a). HTC 
decreases 30 times or more when AEp changes from zero 
to 10T. The latter value is realized for a reasonable differ- 
ence in the carrier density around An = 5 • 10 12 cm -2 at 
T = 300° K. Maximal value of AR varies to some degree 
with the relaxation frequency 7 in graphene as shown in 
the right inset and with the distance d as shown in the 
left inset. It has to be stressed that AT? between iden- 
tical inactive dielectrics covered with graphene is larger 
than that between two Si02 plates. Note that Si02 is 
considered as one of the very best material for the ra- 
diative heat transfer, for which Rsio 2 — 296 W/m 2 K at 
T = 300° K. With temperature decrease Rsi0 2 decreases 
fast due to the thermal factor in ([T]), but AT? decreases 
much slower because of plasmon tuning. As the result the 
ratio AR/ Rsio 2 increases significantly with T decrease 
as shown in Fig. EJb). 

In conclusion, we analyzed the change in the near- 
field radiative heat transfer for materials covered with 
graphene. Plasmon excitations in graphene change dras- 
tically material's ability to the heat transfer. Plasmon 
frequency is tunable and depends on the wavenumber, 
Fermi level, temperature, and the substrate dielectric 
constant. Plasmons can be tuned with the resonances in 
the opposite body to maximize the HTC. The strongest 
effect is reached for terahertz- inactive dielectrics with low 
dielectric constant covered with graphene. In this case 
the HTC is larger than for the best known materials es- 
pecially at low temperatures. HTC can be reduced 100 
times or so by changing the relative carrier concentra- 
tion in different bodies. This opens up the possibility to 
control the heat flux at frequencies as high as 100 GHz. 



[1] D. Polder and M. Van Hove, Phys. Rev. B 4, 3303 (1971). 
[2] J. -P. Mulet, K. Joulain, R. Carminati, and J. -J. Greffet, 

Appl. Phys. Lett. 78, 2931 (2001). 
[3] A. Kittel, W. Muller-Hirsch, J. Parisi, S. A. Biehs, D. 

Reddig, M. Holthaus, Phys. Rev. Lett. 95, 224301 (2005). 
[4] K. Joulain, J. -P. Mulet, F. Marquier, R. Carminati, and 

J.- J. Greffet, Surf. Sci. Rep. 57, 59 (2005); A. I. Volokitin 

and B. N. J. Persson, Rev. Mod. Phys. 79, 1291 (2007). 
[5] S. Shen, A. Narayanaswamy, and G. Chen, Nano Lett. 9, 

2909 (2009); E. Rousseau, A. Siria, G. Jourdan, S. Volz, 

F. Comin, J. Chevrier, and J.-J. Greffet, Nat. Photonics 

3, 514 (2009). 

[6] M. Laroche, R. Carminati, and J.-J. Greffet, J. Appl. 

Phys. 100, 063704 (2006); K. Park, S. Basu, W. P. King, 

and Z. M. Zhang, JQSRT 109, 305 (2008). 
[7] Z. W. Liu, Q. H. Wei, and X. Zhang, Nano Lett. 5, 957 

(2005); L. Wang, S. M. Uppuluri, E. X. Jin, and X. F. 

Xu, Nano Lett. 6, 361 (2006); B. J. Lee, Y. B. Chen, and 

Z. M. Zhang, JQSRT 109, 608 (2008). 
[8] Y. De Wilde et al., Nature 444, 740 (2006). 
[9] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 



(2009). 

[10] K. I. Boloting et al, Solid State Commun. 146, 351 
(2008). 

[11] Y.-M. Lin et al, Science 332, 1294 (2011). 

[12] M. Jablan, H. Buljan,and M. Soljacic, Phys. Rev. B 80, 

245435 (2009); F. H. Koppens, D. E. Chang, and F. J. 

G. de Abajo, arXiv: 1104.2068; L. Ju et al, Nat. Nan- 

otechnol 6, 630 (2011); 
[13] B. N. J. Persson and H. Ueba, J. Phys.: Condens. Matter 

22, 462201 (2010); A. I. Volokitin and B. N. J. Persson, 

Phys. Rev. B 83, 241407(R) (2011). 
[14] B. Wunsch, T. Stauber, and F. Guinea, New. J. Phys. 8, 

318 (2006); E. H. Hwang and S. Das Sarma, Phys. Rev. 

B 75, 205418 (2007). 
[15] O. Vafek, Phys. Rev. Lett. 97, 266406 (2006). 
[16] G. Gomez-Santos, Phys. Rev. B 80, 245424 (2009). 
[17] V.B. Svetovoy, Z. Moktadir, M. C. Elwenspoek, and H. 

Mizuta, EPL 96, 14006 (2011). 
[18] L. A. Falkovsky and S. S. Pershoguba, Phys. Rev. B 76, 

153410 (2007); T. Stauber, N. M. R. Peres, and A. K. 

Geim, Phys. Rev. B 78, 085432 (2008). 
[19] I. V. Fialkovsky, V. N. Marachevsky, and D. V. Vassile- 



5 



vich, Phys. Rev. B 84, 035446 (2011). 
[20] V. N. Kotov, B. Ucho a, V. M. Per eira, A. H. Castro Neto, 

and F. Guinea. larXiv:1012. 34841 
[21] N. D. Mermin, Phys. Rev. B 1, 2362 (1970). 
[22] P. J. van Zwol, K. Joulain, P. Ben Abdallah, J. J. Greffet 



and J. Chevrier, Phys. Rev. B 83, 201404(R) (2011). 
[23] Y.-M. Lin et al, Science 327, 662 (2010); Y. Wu et al., 

Nature 472, 74 (2011). 
[24] X.-Y. Zhao and H.-J. Liu, Polym. Int. 59, 597 (2010). 



